Vibrational Study of 13 C-enriched Ceo Crystals 
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Abstract 

The infrared (IR) spectrum of solid Cgo exhibits many weak vibrational 
modes. Symmetry breaking due to 13 C isotopes provides a possible route 
for optically activating IR-silent vibrational modes. Experimental spectra 



> 

(N 

m 

o 



ctf 



O 



and a semi-empirical theory on natural abundance and J C-enriched single 



'-£} ■ crystals of Cgo are presented. By comparing the experimental results with 

the theoretical results, we exclude this isotopic activation mechanism from 



the explanation for weakly active fundamentals in the spectra. 
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The availability of high quality Cgo single crystals has provided an opportunity for the 
detailed study of weak vibrational features in the infrared (IR) and Raman spectra of this 
molecular solid. Recent studies identified some of these modes as being weakly activated 
symmetry-forbidden modes in IRUta and Ramanffi measurements. Second-order combina- 
tion modes were also observed in these spectra. These results enabled the authors to use 
group theoretical arguments to extract candidates for all 46 vibrational modes of Cgo from 
the experimental dataBs 

Since the icosahedral (Ih) symmetry of the isolated fullerene molecule allows for only 
four IR active modes, this symmetry must be broken by some agent. Possible mechanisms 
for activating such "silent" modes include 13 C isotopic impurities, Van der Waals forces 
from other Cgo molecules (crystal field effects), anharmonicity, electric field gradients at 
surface boundaries, impurities, and dislocations. Several previous worksa'Erera have argued 
that 13 C isotope substitution is a likely agent for activating silent modes since for the natural 
abundance of 13 C, about one half of all Cqq molecules would have at least one 13 C which 
would break the Ih symmetry and activate all the vibrational modes. In the present study, 
we experimentally and theoretically explore this idea by studying Ceo crystals made from 
natural abundance and 13 C-enriched C 60 . 

The natural abundance 13 C single crystals were sublimated as described previously^ 
from well cleaned Cgo powder commercially produced.0 Typical crystals have dimensions 
~ 1 x 1 x 0.5mm. The 8% 13 C-enriched crystals were sublimated from enriched soot pre- 
pared as reported earlierQ and their chemical and crystallographic quality were checked by 
high-resolution 13 C solid state nuclear magnetic resonance (NMR)H The isotope enrichment 
was determined by mass spectrometry, showing that the mass distribution follows a Pois- 
son function. A Bomem MB-155 FTIR spectrometer was used to obtain 2cm -1 resolution 
IR spectra in the 400-5000cm _1 energy range. The infrared spectra presented here were 
obtained at room temperature with the crystals mounted over a small aperture with a min- 
imal amount of silver paste. The spectra are referenced to spectra obtained through the 
empty aperture. The IR spectra of such thick crystals has been shown to be independent of 
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atmospherea, so for experimental simplicity the crystals in this study were exposed to air. 

The probability that a Ceo molecule has n 13 C atoms is given by the Poisson distributions, 
p n , in agreement with the mass spectrometry of the 8%-enriched 13 C samples.^ This function 
is presented in Figure [I] for Ceo molecules made from natural abundance carbon (1.108% 
13 C) and from 8% 13 C-enriched carbon. The natural abundance probability curve shows 
that about one third of the molecules will have exactly one 13 C atom which will weakly 
break the icosohedral symmetry of pure 12 Ceo- From the 8% probability function we learn 
that only about 1% of Ceo molecules will be pure 12 Ceo; most will have 3 to 7 13 C atoms per 
fullerene. 

The IR data are presented in Figure |2[ The uppermost curve is the transmission of 
a natural abundance C 6 o single crystal, as published previouslycl The four IR-active F lu 
vibrational modes are observed to be completely saturated around their known positions 
of 527, 576, 1182, and 1427cm -1 . All other observed absorptions are due to IR-forbidden 
fundamental vibrations or second (or higher) order modesHfil™ which have become weakly 
active in the crystal. One of the weak modes at 710cm- was fitted in a previous study!; its 
intensity (proportional to the square of the plasma frequency of the Lorentz oscillator, uJ^) 
is approximately one order of magnitude smaller than that of a typical F lu mode. 

The middle curve of Figure ||] displays the IR spectrum of an 8% 13 C-enriched Ceo crystal. 
When comparing this spectrum with that of the natural abundance crystal (top curve of Fig. 
0), the primary changes observed are that all of the visible modes have broadened and slightly 
softened. Note, however, that the relative intensities of the weakly active fundamental and 
higher order vibrational modes do not change significantly. Nearly all the features of the 8% 
enriched sample transmission spectrum can be modeled by simply broadening the spectrum 
obtained for the natural abundance crystal using the functions from Figure [I] and assuming 
the energy of a mode is proportional to 1/ ' \fm (like in a simple "ball and spring" model). 
(This relation between mode broadening and mass distribution was described in greater 
detail for the ^4 S (2) Raman fundamental by Guha et a/..M) Carrying this procedure out by 
applying the 8% Poisson distribution point by point to the top spectra of Fig. ||, we obtain 
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the bottom curve of Figure |2]. A careful comparison of the measured 8% enriched spectrum 
and the spectrum generated from the broadening reveals that they are indeed very well 
matched, indicating that 13 C substitutions change the frequency of the mode, but do not 
influence the intensity. Raman spectra obtained on these samples in which forbidden and 
higher order modes are observed give the same results.EJ 

If the fundamental modes are being activated by symmetry breaking due to isotopic 
substitutions, one would expect the 13 C-enriched samples to have greater symmetry break- 
ing, and therefore greater IR-activity for originally silent modes. To test this, we used a 
theoretical model of Cgo based on the bond-charge model of Onida and Benedekfca (the only 
change was to use the better experimental valuest3 of 1.4A for the hexagon-hexagon (6-6) 
bond length and 1.45A for the pentagon-hexagon (5-6) bond length). The vibrational fre- 
quencies calculated with this model are within 3% of the experimental IR and Raman mode 
frequencies. 

The calculation of the IR intensity of a given mode requires, in principle, the knowledge 
of the electronic eigenstates of the deformed molecule. To avoid the numerical complexity of 
such an approach we used an approximate procedure, based on the experimentally known IR 
intensities of the four allowed F\ u modes and on the calculated vibrational eigenvectors of the 
isotopically changed molecule. When expanded in terms of the original eigenvectors of the 
pure 12 C molecule, some or all of the new eigenvectors had components in the directions of 
the original F\ u eigenvectors. Each degenerate species of the unperturbed F\ u eigenvectors 
are chosen to carry their dipole moments in the Cartesian x, y, and z directions. The 
intensity of each of the four IR-active F\ u modes was then taken from experiments and are 
denoted A±, A 2) A3, and A*. The IR-intensity of mode I is then calculated using 

where Q}/ is the component of eigenmode I in the direction of the the i th degenerate mode 
of the original F\ u (j) eigenvector; i sums over the triple degeneracies and j sums over the 
four Fi u modes. The isotopic content is modeled by a Monte Carlo simulation: each of the 
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60 carbon atoms is given a 13 C or a 12 C mass with appropriate probabilities representing the 
natural abundance or 13 C-enrichment (1.108% and 8% 13 C, respectively). The dynamical 
matrix for the resultant molecule was then diagonalized to calculate the eigenvectors, and 
Eqn. [J] was applied. For a given isotope concentration, this process is repeated until a 
general sampling of all possible 13 C configurations is calculated and the results are averaged 
together. For the natural abundance case, there are relatively few possible configurations, 
but for the 8% 13 C-enriched case, there are a huge number of configurations. We sampled 
over 15,000 possible configurations for the 8% case, which should suffice to obtain the correct 
trends. 

We plot the resultant calculated IR- intensity for 1.108% and 8% 13 C Ceo molecules in 
Figure [3[ The Lorentzian widths of all modes have been taken to be uniformly 3cm -1 . The 
four Fi u modes in Figure |3] have peak heights of approximately 2000 {u^ ~ 5000cm -1 ), 
so the intensities of the weakly activated modes are approximately a factor of 10 4 smaller. 
This is in contrast to the experimentally observed intensity ratio of only a factor of ~ 10. 
Furthermore, the average intensity of the weakly activated modes increases by a factor of 
4.45 in the 8% enriched calculation compared to the natural abundance calculation, but no 
similar increase was seen in the experiment. Very similar theoretical results are given in a 
recent detailed calculation^ of the effects of isotopes on the vibrations of Cqq. 

We must therefore conclude that isotopic activation of previously silent modes is not the 
symmetry-breaking mechanism in Cgo- 

Since the isotope effect is eliminated, the forbidden modes must be activated by another 
symmetry breaking mechanism. Surface effects can be ruled out since the phenomenon is 
independent of sample thickness. Defects in the bulk, like dislocations or vacancies are 
also less likely, since the measurements on samples under pressure, or samples subjected 
to plastic deformation by pressure do not exhibit a dramatic change in these weak linesB 
Oxygen bound to the fullerenes has been excluded since the number and strength of the 
weak vibrational modes are unchanged when obtained from Cgo crystals that have never 
been exposed to oxygenB Impurities in the crystals are a possibility, however we point out 
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that crystals made from different C 60 sources and using various growing conditions result 
in exactly identical IR spectra in every feature.ETH A remaining serious candidate is the 
symmetry breaking due to the first neighbor interactions between the C^o molecules in the 
solid. 

The first question is, how is it possible that most of the weak fundamentals are visible 
at room temperature, when the X-ray structure is face centered cubic, with one fullerene 
molecule per unit cell.EJ In this state the molecules are rotating rapidly and the symmetry 
breaking crystal fields fluctuate with this rotation. To clarify this issue, we performed a 
simple model calculation with a fluctuating external potential breaking the symmetry of an 
otherwise symmetric arrangement of atomsHI All other interactions with the surroundings 
were represented by a damping constant. We found that a forbidden mode is activated 
{i.e. the system absorbs energy from the external oscillating electric field in a resonance like 
fashion) as long as the correlation time of the fluctuations is not shorter than the lifetime 
of the resonance (which, in our calculation, was determined by the damping constant). 

To apply this idea to the present case, we need to know the natural lifetime of a par- 
ticular vibrational mode, and the correlation time of the fluctuations of neighbors. The 
first parameter can be estimated from the linewidth of the well resolved F lu vibrational 
modes as measured in detail by Homes et a/.Jia To avoid complications due to the "inho- 
mogeneous broadening" we take a linewidth from the measurement below the orientational 
ordering transition temperature (7 = 0.5cm -1 ) and obtain a lifetime of 60psec. At higher 
temperatures the lifetime decreases, therefore this is an upper limit. On the other hand 
the rotational diffusion constant was determined by neutron scattering measurements as 
Dr = 1.4 x 10 10 sec _1 at T=260K, right above the rotational phase transition.!^ This implies 
that a 60 degree rotation takes about 130psec, which is indeed longer than the upper limit 
to our estimated resonance lifetime. 

According the the neutron scattering data, the rotational diffusion constant increases by 
a factor of two, to 2.8 x 10 10 sec _1 at 520K. However, the lifetime of the resonance is also 
expected to decrease. In an insulator, the primary channel of phonon decay is due to the 
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interaction between the phonons, and the enhanced population of the phonon states leads 
to shorter lifetimes at higher temperatures. These considerations lead to an experimentally 
verifiable prediction: at higher temperatures the "forbidden" lines should either disappear 
(if the lifetime of the resonance has a weak temperature dependence) or the integrated 
intensities will not change while the lines broaden (so that the lifetime stays shorter than 
the fluctuation time). To test this prediction, we measured the IR spectrum of a crystal at 
elevated temperatures. All observed modes do broaden without a significant change in line 
intensities up to a temperature of 650K. 

Below T~ 250K, the orientational ordering transition temperature, the halting of free 
rotation of the Ceo molecules leads to a crystal structure with four fullerenes in the unit 
cell.E9 In this structure all the ungerade modes of the molecule are allowed in the IR by 
symmetryO. The ideal Pa3 structure retains a center of inversion symmetry at the center 
of each Ceo molecule and thus the gerade modes are still IR-forbidden. However there is a 
"defect orientation" in the actual structure which will break this inversion symmetry in a 
significant portion of the Ceo molecules.^ 

According to these arguments there is reason to believe that the nearest neighbor sur- 
roundings of the molecules provide the symmetry breaking necessary for the activation of 
the forbidden modes. This mechanism is supported by the fact that many of the weakly 
active modes exhibit a sharp change when cooled through the orientational ordering phase 
transition.BH On the other hand the characteristic temperature dependence of the "defect 
orientation" fractionE3 should lead to a gradual decrease in the intensity of the weakly active 
IR lines originating from the gerade modes (in particular Raman modes), and this was not 
seen in our experiments.0 

To summarize, we have shown that isotopic symmetry reduction is not the activating 
mechanism for symmetry forbidden vibrational modes in Ceo as has been previously pro- 
posed. We have shown that crystal field effects can explain the weak modes both above 
and below the orientational ordering transition temperature, but further experimental and 
theoretical work is required to fully understand the nature of the symmetry breaking in 
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fullerene crystals. 
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FIGURES 
FIG. 1. Plot of the probability function that a Cgo molecule has n 13 C and (60-n) 12 C atoms 
for the natural abundance of 13 C and for 8% enriched samples as measured in the present study. 

FIG. 2. IR spectra of Cgo single crystals; the two upper curves have been vertically offset for 
clarity. The top spectra is from a crystal made from natural abundance 13 C Cgo- The middle 
curve is the spectra of an 8% 13 C-enriched Cgo crystal. The bottom curve is generated from the 
top curve using a simple model of frequency shift and the mass distribution given from Figure 1 
for 8% enrichment. 

FIG. 3. Calculated IR intensities of vibrational modes in Cqq with 1.108% and 8% isotopic 
substitution. The four IR-active F\ u modes have intensities of ~ 2000 which are well off scale; the 
scale is chosen so that the weakly activated symmetry-forbidden modes can be seen. 



